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1.  INTRODUCTION 


This  report  examines  Gough’s  (1974)  derivation  of  the  solid-phase  portion  of  the  Navier-Stokes 
equations  of  motion  for  the  NOVA  series  of  interior  ballistic  (IB)  codes.  The  motivation  for  this  report 
is  the  better  understanding  of  the  assumptions  used  in  the  formulation  and  how  these  may  pertain  to 
experimental  data  input  for  constitutive  modeling.  Also,  this  report  provides  the  derivation  in  a  concise 
location  and  step  by  step  format  for  the  unfamiliar  reader. 

2.  DERIVATION  OF  THE  CURRENT  SOLID-PHASE  EQUATIONS 

The  following  is  the  derivation  of  the  solid-phase  Navier-Stokes  equations  for  IB  using  the  original 
formal  averaging  procedure  used  by  Gough  in  1974.  Gough  did  not  explicitly  demonstrate  this  derivation 
of  the  solid-phase  portion  in  his  thesis.  The  more  general,  one-dimensional,  isothermal  Navier-Stokes 
equations  include  the  solid-phase  continuity  equation. 

^  V-(p^,)  =  (1) 


and  solid-phase  momentum  equation, 

9(p^  ) 

_^.v.(p,«^«^.r,)=4»,v,.  (2) 

The  principle  assumptions  to  be  made  initially  are:  (1)  that  the  propellant  is  compacted  isothermally: 
(2)  the  stress  tensor  is  macroscopically  isotropic;  and  (3)  the  .solid-phase  mass  production  term  is  zero. 
The  isothermal  assumption  considering  the  thermal  scale,  including  combustion,  assumes  that  the 
temperature  rise  of  the  propellant  due  to  compaction  is  small  with  respect  to  the  energy  released  by 
combustion.  The  isotropic  assumption  was  made  so  the  mixture  stress  (i.e.,  the  intrinsic  stress  multiplied 
by  the  solid-phase  volume  fraction)  could  be  more  easily  isolated  from  the  gas  pressure;  this  also 
eliminates  any  shear.  The  mass  production  is  assumed  zero  over  small  time  intervals  due  to  the  small 
production  losses  in  comparison  to  the  convective  term. 

Following  the  formal  averaging  procedure  outlined  in  Gough’s  thesis,  the  momentum  equation  may 
be  reduced  through  the  following  steps.  No  solid-phase  mass  production  implies 


1 


(3) 


3(p,«j)  „  _ 

— gj —  +  +  r,)  »  0. 


Turbulent  fluctuations  of  a  variable  arc  accounted  for  through  the  following  typical  notation,  i.e., 


<9y  *  p* 

(4) 

<u  u* 

(5) 

and  so  on  for  other  variables.  Formal  averaging  may  be  written  in  either  intrinsic  (porosity  weighted)  or 
iK)nintrinsic  forms. 


<pj(x,t)u  ^(X,t)u  *  T^(x.t)> 

»  ejix,t)<p^(x,t)u j(x,t)u j(x,{)  +  T^(x,t)>^ 
or  =  r  g(y"  -  x,x  -  i)[p^(y,x)Uj(y,x)u,(y,x)  +  T ^{y,x))dVdx  .  (6) 

S  S 


The  intrinsic  average  is  identified  by  the  term  with  the  superscript  5.  Equation  6  introduces  an  arbitrary 
weighting  function,  g.  which  is  assumed  to  be  continuously  differentiable  and  satisfies  the  normalization 
condition 


JJJ/ 


g{y,x)dVdx  «  I 


(7) 


Using 


J  -  x,x  -  t)) 

5^) 


3(g0'  -  x,x  -  t)) 


(8) 


the  development  of  the  divergence  term  in  Equation  3  may  commence  with  the  knowledge  that  p,  f,  and 
u  are  independent  of  the  spacial  dummy  variable,  y,,  allowing 
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d 

dx: 


<pjix,t)u  ^(x,t)u  j(x,t)  +  T^(x,t)> 


Lx  s(y  -  ~  0-3i(pj(y,'i)«i(y.x)M  ,(/ ,t)  +  f  jiy,x))dVd  X 

~  j^x  -  X.x  -  t)(p^(y.x)u^(y.x)u,(y.x)  *  f,(y.x)j 


dVdx  .  (9) 


Using  the  definition  of  the  nonintrinsic  average  and  Gauss’  theorem,  both  the  divergence  and  temporal 
derivative  tenns  of  the  momentum  equation  may  be  expanded  as; 


^^<p^ix,t)u^ix,t)u  /,x,t)  +  T^(x.t)> 

<^[p^(x,t)u^(x,t)u^(x,t)  *■  f^(x,t))> 
ax, 

-  -  X,x  -  t){p^(y,x)u^(y,x)u^(y,x)  ->■  f  ,(y.x))]  ’  ndA 


(10) 


and 


a 

■5? 


<p,(x,r)M^(x,r)> 


+ 


<  -g-  |pj(x,r)M^(x,r)j> 
Op,(y,x)«j(y.'c))  w  •  ndA  , 


(II) 


where  w  is  the  relative  velocity  between  the  phases.  The  following  nonintrinsic  averaged  form  of  the 
momentum  equation  may  then  be  written  using  Equations  10  and  11, 


^<P,«f> 


’'S 

TT 


■{<p,«,« 


s> 


.  ( M  - w ) ) • ndA 


(12) 
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Allying  ihe  tnaimic  average  to  ^iuMion  12  ^oduces 


•  V.  U(«p, 


^  ( r ,  *  P,  “ ,  <  “ ,  "  ♦*■  )  >  * 


<J3i 


To  Auther  nech^  the  terms  of  tius  ettimkm.  one  may  wnte  out  the  etpwuuons  for  the  svmf»| 
tecfmique  fm'  products  of  two  vanabkes. 


<  P,  •  <(<  P,  >'  •  P,*)  (<  >'  •  «,*)>' 

J  I  **  1  •  *  #  #  *•  •  \  ^  n4  i 

m  <1<  p^>  <tf,>  *  •  <  p^>  il  ^  /  j  ^  *  ' 


and  three  variables. 


/  s  *  \ 

\<Px^  *  P,  j 

P“P  "“xl 

|<«4,>‘  ♦  «,  /> 

(S  ^  i  •"*  f  t  “’♦  I 

<p^><«,>  ♦  <P,  >«,  *P,  <“,>  *P,  *  “»P 

•  <(<p^>*<«^>^<M^>'  ♦  <P^>''«*<«^>’  *  p*  <  14  ^  >'<  >'  •  P*tf* 

♦  <p,  >  <«,>«,  *  <P^>  ’  U  ^  ♦P^  •  p#“4“* 


<  «  ^  >' 


(55) 


The  average  of  a  ntK:tuating  term  alone,  or  the  average  of  a  fluctuating  term  multiplied  by  an  averaged 
term,  is  zero,  allowing  the  terms  containing  a  single,  averaged  fluctuating  variable  to  be  negleaol.  The 
terms  ctmtaining  two  or  nrore  fluctuating  variables  multifriied  together  and  then  averaged  are  not  zero. 
It  is  also  notol  that  die  avera^  of  an  averaged  properry  is  the  average,  tlm  is. 
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<  <  p,  «,>  *>  -  <  P,  Ut>* 


( lO) 


This  resulis  in  averageU  mocfnentum  serais  siniilaf  lo  Gougti  s  equaitsiri  i3  1  4  6|. 


<  p, « ,  > '  -  <  p,  > '  <  M  ,  > ' 


J 

> 


(IT) 


and  equation  (3.i.4  7), 


<P,  U,M^>  «<pj>  <«,>  <»»,>  •<P,  > 


»2<«,>  <p,  u,>  •<P,  • 


for  the  gas  phase  terms. 

With  this  genera]  background,  the  averaging  of  the  mornenturn  equation  may  proceed  The  temporaj 
term  is  expanded  using  Equation  17  as 

-  a(e,<p,>^)  a(e,< 

- , -  -  <H>  - - -  »  -  *  . --r - 

df  dt  dt  ^ 

The  divcrgctKc  term  of  Equation  12  expands  to 


(M) 


V  •  I  <  pj  y  <  u ,  >*  <  M ,  ♦  <  p,  >^  <  « , « ,  >* 

2<«,>"  <p,  «,>  >  <p,  «,«,>  I 
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Using  the  relMion 

V  ‘  {u  It)  •  (V  '  u)  u  *  It  ■  (V  u) .  (21) 

which  may  be  easily  (kmonstnted  in  Canesi»J  coordinates  for  a  wed  bdwved  vmcNr  m .  the  csssmston 
of  the  divergence  term  given  in  E^taiion  20  may  (hen  be  wnoen  similar  id  Googh's  oquautm  <3.2.3  4) 
for  the  gas  phase: 

♦  <«^>*  V  •  *  e,<P*,“*,>^) 

♦  e,  <  p,  >^  <«  ,>^  •  V  <  a ,  >'  *  V 

♦  V  •  (€,(<p,>^<tr*,«*^>'  •  <a,>'<P*,«*,>'  *  <P“,*‘*,“*,>lj  (22) 

With  this  and  the  temporal  term  in  Et|uation  19  substituted  into  Ecfuation  12.  the  mocnetttum  equation 
becomes 

<“> - 57 -  ‘  --57”  * - 57 - 

♦  e,<p*^«*^>^  ■  V<«^>'  ♦  <«^>*  V  •  |c,<p,>*<ir^>*  •  t,<P*, “*,>*) 

♦  V  ■  (c^(<p^>^<tt*,«*^>^  ♦  <i7^>''<p*,«*^>^  ♦  <P* jfH* ,>*)] 

•  *  P,U,(U,  -  W))-  ndA  (23) 
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Suniiaf  opemions  with  the  fcsthctions  of  constant  density,  and  nomeacung  interfaces,  transform  the 
cotuinuity  equation  (Equation  1)  into 


d(e,<p,>") 


♦  V  •£,(  <p^>^<M^>*)  •  0  . 


(24) 


Recognizing 


D<u  >*  h<u>*  I  I 

c,  <p,>^ — f —  -  e,<p,>" — —  ♦  C,<p,>"(<«,>'-V)<«,>' 


(25) 


and  a{^lying  Equation  24.  with  the  assumption  of  no  relative  motion  of  the  p^tases  (i.e.,  u,  »  Ug).  to 
Equation  23  results  in 


D<u  ^>‘ 
_ 


(26) 


where 


e 


2x  • 


(27) 


Transformation  of  the  unknown  solid  phase  surface  integral  into  a  gas  phase  boimdary  condition  may  be 
accomplished  by  introducing  an  intrinsic  intergranular  stress  describing  the  total  stress  on  the  propellant 
bed  as 


<f,>^MR*<fg>»,  (28) 

where  the  gas  phase  contribution  is  considered  as  a  hydrostatic  pressure  only 

<f,>*  -  •/  .  (29) 

and  the  solid  gas  phase  stress  dinerence  as 

if,-fg)  ■  n  -nap,  (30) 
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where  this  interphase  stress  difference  is  usually  zero  unless  the  system  is  untter  reaction  or  surface 
tension.  Using  the  following  integral  theorem  for  the  averaged  value  of  aas  intrinsic  property, 


J  g  dVdt  -  -V(e,<f^  . 


(31) 


ss 


and  applying  Gauss’  theorem  to  this  provides 


g<T,^ 


ndA  »  -  <  r j  5^  •  VCjj  . 


(32) 


thus  allowing  the  integral  term  from  Equation  26  to  be  written  using  Equations  30  and  32  as 


<f,>»  -Ve,  .  j^^^g  (f,  -  <f,>>)-ndA 
*  (&P)’  ndA  . 


(33) 


The  coupling  stress  integral  on  the  left-hand  side  of  Equation  33  is  of  both  the  averaged  and  fluctuating 
portion  of  the  solid  phase  stress  tensor.  The  first  integral  on  the  right-hand  side  of  Equation  33  is  merely 
the  integral  of  the  fluctuating  portion  of  the  gas  phase  stress  tensor  given  by  the  relationship 

Equation  26  may  be  rewritten  using  Equation  I  and  the  integral  form  of  Equation  32,  noting  that  tl^re 
is  a  sign  convention  reversal  of  n  for  agreement  with  the  gas  phase,  as 

e,<p,>-£^.V.(e.<f,>*) 

»  02,  -  <Tg>^  -Ve^  +  4)*«^^**  ~ 

Using  Equation  28  to  expand  the  ^ss  tenn  and  letting 

<Ps>'^Pp’  (36) 
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results  in 


^  D<u  /  “  1 

-  9l,  >  ■<■<«-  4A.9e,  ,  ,37) 

where  the  subscript  p  denwcs  a  property  of  the  particle  The  panicle  imcrfacitl  force  term  given  by  the 
integral  ot  the  right-hand  side  of  Equation  37  is  itprcscracd  as 


■  ndA  » 


<F>'  . 


(3«) 


where  <F>‘  is  the  microscopic  intcraaion.  drag,  between  the  media.  Suhaituting  this  into  Equation  35 
results  in  the  following  form  of  the  solid  f^ase  momentum  equation: 


^s'^Pg  *  ^  "/  • 


(39) 


where 


D  d 

^  **5?: 

p  p 


♦  u 


V 


(40) 


and  the  inteq^ase  forces,  as  well  as  tlK  turbulent  momentum  source  term,  are  combined  as 


/ 


02»  "■ 


hh^<F>'  -  ApVe, 


(41) 


Further  manipulation  of  the  averaged  equations  reformulates  the  momentum  equation,  after  linearization 
of  the  stress  term,  into  a  recognizable  wave  equatitm  as  follows.  The  following  averaged  form  of  the  solid 
phase  continuity  equation. 
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(42) 


de 


P/rgy^  ^ 


may  be  rewritten  as 


r,  - 

■  P;>^x^'“p  •  • 

V 


(43) 


where 


^  -  4-  ^  a;  *  V  • 

Ot^  IT 


(44) 


Expanding  the  stress  term  in  the  momentum  Equation  39,  with  R  now  represented  as  a  hydrostatic  stress. 


V  ( e,/? )  »  e,  V  /?  +  /?  V  e,  . 


(45) 


and  further  expanding  the  gradient  term. 


(46) 


where  the  mixture  stress  of  the  particles  is  o  =  e^,  provides 


V(op 


(47) 


Let  G(eJ  be  a  compressive  modulus  such  that 


V(o  )  -  C?(e,)Ve 


(48) 
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Substituting  this  into  Equation  39  results  in  the  following  nonlinear  momentum  equation  to  be  solved  in 
the  IB  code 

+  <3(e,)Ve^  -/  .  (49) 

P 

Further  assumptions  are  now  made  in  order  to  model  the  stress  term  in  the  momentum  equation.  These 
assum^Mions  consist  of;  [4]/=  0,  no  fluids  in  the  propellant  bed;  [5]  «  •  v  =  0,  the  particles  are  initially 
at  rest  (i.e.,  for  small  perturbations  about  a  quiescent  state);  and  [6]  /?  =  /?(ep  <mly.  The  assumfxion  of 
no  fluid  in  the  bed  allows  the  gas  pressure  term  to  be  i.cglected  when  computing  the  mixture  stress.  The 
small  perturbations  assumption  eliminates  the  convective  portion  of  the  substantial  derivative,  thus  the 
particles  do  not  "flow,"  thus  only  accelerations  and  interaction  forces  exist.  Assumption  [6]  is  only  for 
simplicity;  of  course,  1?  is  a  function  of  temperature,  contact  surface  area,  and  many  otter  things  which 
may  or  may  not  be  experimentally  measurable.  The  one  variable  that  is  measurable  in  a  bed  compaction 
test  is  the  change  in  volume.  This  allows  the  computation  of  the  porosity,  assuming  incompressibility  or 
given  the  intrinsic  equation  of  state  of  the  particles  from  a  second  experiment. 

A(^lying  these  assumptions  reduces  the  momentum  Equation  49  to 


n 


(50) 


and  the  continuity  equation  (Equation  1)  to 


r,  - 

IT  “ 


(51) 


One  approach  to  deriving  the  solid  phase  wave  speed  is  to  take  the  derivative  of  the  momentum  equation 
with  respea  to  time.  The  assumption  that  the  intergranular  stress  /?  is  not  a  function  of  time  permits. 


Pp~ 


dTu 

t 


'SO 


)^Ve,.0 


(32) 
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where 


(53) 


Using  continuity  equation  (Equation  SI),  Equation  S3  may  be  written  as 


(54) 


which  is 


+  Vcuil(curi«^)) 


(55) 


Inotational  flow  implies  that  the  curl  of  the  velocity  is  zero.  Substituting  this  into  the  momentum 
equation,  and  dividing  through  by  the  porosity  and  density  results  in 


v2-  0 

(>P 


(56) 


This  is  a  wave  equation  having  a  velocity  of 


(57) 


where,  for  the  linearized  case. 


Go(£,o) 


(58) 


In  the  general  nonlinear  case,  the  function  G  behaves  as  a  "stifiness"  modulus  which  is  the  product 
of  the  density  and  the  propagation  speed 


Git,) 


t'p  p 


(59) 
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3.  RHEOLOGICAL  REQUIREMENTS 


In  the  theory  underlying  Gough's  XNOVAKTC  code  (1990),  the  constitutive  assumption  which 
defines  the  solid-phase  stress  tensor  is  embedded  into  the  function  G(e^)  by  specifying  a  functional 
dependence  of  the  propagation  speed  which  is  assigned  for  a  system  undergoing  loading  as 


(60) 


and  unloading 


(61) 


where  the  user-supplied  constant  Uj  represents  the  speed  of  propagation  during  compressive  loading  when 
the  bed  is  at  the  settling  porosity,  and  the  consumt  02  represents  the  propagation  speed  derived  from 
the  modulus  during  unloading  from  any  state. 


A  functional  dependance  of  aJe.)  may  be  developed  by  rewriting  Equation  50  as 

r  .  ^2 


de. 


=  G(eJ  =  - 


80 


) 


de„ 


(62) 


Integrating  this  relationship  from  to  where  is  the  current  gas  phase  porosity  and  is  the 
settling  porosity,  and  a^  is  the  user-input  propagation  velocity  of  an  infinitesimally  small  disturbance, 
generates  the  loading  function  used  within  the  XNOVAKTC  IB  code. 


o/ep 


g  go 
60 


(63) 


Both  Oi  and  02  are  experimentally  determined  constants.  The  speed  of  propagation  of  a  small  disturbance 
up  is  usually  about  an  order  of  magnitude  smaller  than  that  of  02^.  These  relations  can  be  adjusted  to 
the  Robbins  and  Conroy  (1991)  model,  which  was  performed  on  the  data  from  the  Birkett  (1981)  test  rig 
as  shown  in  Figure  1. 


The  four  experimental  data  sets  in  Figure  1  end  at  different  final  force  levels,  while  the  functional  data 
are  plotted  for  four  different  disturbance  propagation  speeds,  aj,  increasing  from  100  m/s  to  4(X)  m/s  from 
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the  lower  to  the  upper  computed  curve,  is  0.45  and  the  diameter  of  the  test  cylinder  was  7.62  cm. 
The  curves  shown  in  Figure  1  exhibit  the  correct  trend.  However,  the  chosen  fixed  wave  speed  does  not 
represent  the  curving  of  the  propellant  experimental  data  as  well  as  one  might  want. 

TTie  experiments  which  determine  the  user-input  propagation  velocity  usually  choose  one  value  of  the 
modulus  from  limited,  typically  below  1  MPa,  axial  stress  levels.  The  data  are  obtained  by  compressing 
a  bed  of  propellant  in  a  cylinder  and  measuring  the  applied  force  and  corresponding  displacement.  From 
this,  a  modulus,  G,  is  determined  at  the  settling  porosity  and  then  the  simple  relation 

fli  - 

is  used  to  obtain  the  compressional  propagation  speed  for  a  small  disturbance.  Currently,  only  a  limited 
amount  of  data  is  available  for  gun  propellants  of  interest  at  this  time  (Conroy  1992).  In  addition, 
estimating  the  value  of  the  modulus  of  the  propellant  involves  a  certain  amount  of  ambiguity,  because  the 
0.2%  rule  for  determining  the  modulus  is  not  applicable  to  a  structure,  such  as  a  bed.  Unfortunately,  the 
predicted  values  of  stress  (particularly  at  small  values  of  porosity,  ep  are  rather  sensitive  to  the  value  of 
ai  (as  shown  in  Figure  1).  which  is  dependent  on  the  calculation  of  the  modulus. 

Class  3  HMX  is  a  granular  explosive  with  a  settling  porosity  of  34.7%.  whereas  TS-3659  is  a  double¬ 
base  ball  propellant  with  a  settling  porosity  of  40.0%.  Figure  2  shows  the  region  in  which  the  propellant 
exhibits  a  stiffening  behavior  as  the  porosity  decreases.  Thus,  to  model  the  entire  range  of  the  propellant 
behavior,  as  potentially  predicted  in  Figure  3  by  the  XNOVAKTC  code,  the  experimental  data  should 
exceed  the  "knee"  effect  which  typically  occurs  at  a  mixture  stress  (bed  stress)  of  about  30  to  40  MPa  for 
double-base  propellants. 

As  a  means  of  improving  the  current  formulation,  at  least  two  approaches  may  be  appropriate.  The 
first,  a  direa  correlation  of  experimental  data,  can  produce  the  loading  function  for  insertion  into  the  code 
and  thus  a  wave  speed  from  the  derivative  thereof.  Second,  modeling  of  the  bed  aggregate  from  single 
grain  behavior  is  more  fimdamental.  In  the  past,  modeling  has  been  performed  on  various  porous 
materials  such  as  the  elaborate  sf^erical  particle  in  the  contact  model  of  Brandt  (1955)  and  the  spherical 
pore  coUapse  model  (Carroll  and  Holt  1972).  Brandt  extends  his  work  to  include  nonspherical  panicles 
continuing  to  apply  Hertz’s  (1881)  theory,  which  is  valid  as  long  as  the  particles  have  an  average  radius 
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70000 


Figure  1.  A  comparison  of  axial  load  predictions  given  bv  Gough  (equation  61)  and  the 

experimental  data  obtained  at  the  Naval  Ordnance  Station  (Biricett  198]~)  for  M30 
propellant  at  219  K  in  a  3-in  diameter  compaction  cylinder. 


Mixture  Density  (g/cc) 

Figure  2.  Experimental  data  for  class  3  HMX  and  TS-3659  and  Kooker’s  correlations. 

Thtfottowing  p*g*  iiutmtiomattf  kft  bkmk. 
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Intentionally  left  blank. 
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of  curvature.  There  is  some  work  on-going  in  this  field  foi  granular  propeUants  at  iht  ARJL  (Ga^onas 
1992).  but  is  tKH  otherwise  being  pursued  to  this  auihor’s  knowledge.  Cosianiirw  ^1983)  gave  a  review 
of  the  various  models  for  energetic  materials  derived  from  ttxk  mechanics. 

A  constimtivc  modeling  comparison  has  been  performed  by  Conroy  and  Kookcr  (1991).  This  report 
highlig^  Kooker’s  (Kooker  and  Anderson  1985)  iiKorporation  of  the  solid  pNhase  equaiitm  of  Hate  to 
model  the  compressibility  and  his  correlation  of  experimental  dau,  a  modified  from  of  Walton's  (1977) 
ecpiation.  This  was  to  demonstrate  the  potential  prediction  (hffetcnces  between  the  current  model  imd  one 
in  which  the  dau  is  correlated  directly,  currently  isotherm  ally.  The  refills  show  some  significaru 
diffeietrces  as  Kooker's  model  tracks  the  solid  material  from  the  sealing  porosity  to  where  it  is  a 
compressed  solid  with  no  voids. 

4.  DISCUSSION  AND  CONCLUSIONS 

The  formulation  of  the  solid-phase  portion  of  the  current  NOVA  family  of  codes  has  been  thoroughly 
described.  The  required  closure  for  this  portion  of  the  IB  model  has  been  exfMained.  Both  analytical 
models,  and  analyiical/cxpcrimcntal  models  have  been  examined  by  researchers  as  a  means  to  provide  the 
stress  strain  relationship  needed  by  the  derived  solid-phase  momentum  equation.  Previous  investigators 
have  resorted  to  using  experimental  propellant  data  to  provide  inputs  for  specific  derived  constitutive 
relationships  from  soil  mechanics.  The  use  of  expciimcnial  dau  directly  has  also  been  investigated  and 
is  a  viable  and  appropriate  surfing  point  for  improving  the  current  model,  however,  the  dau  must  be 
obtained  in  the  regions  to  be  investigated,  which  most  probably  will  include  the  region  above  100  MPa 
as  presented  in  Figure  2.  Noncoincidcntal  efforts  are  currently  underway  to  obuin  quasi-sufic  isothermal 
bed  compaction  dau  to  150  MPa  for  various  propellants  between  BRL  and  Naval  Surface  Warfare  Center, 
White  Oak. 
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LIST  OF  SYMBOLS 


A|  =  User  input  small  distuibancx  propagation  velocity 
a  -  Compressional  wave  speed  used  in  the  solid  phase  wave  equation 
e^  =  Solid  phase  internal  ettergy 

F  =  Microscopic  interaction  bttween  the  media 

f  =  InteT]4iase  forces 

=  Steady-state  drag  per  unit  volume 
G  =  Bed  modulus 

Gg  =  Bed  modulus  at  the  settling  porosity 
g  =  Normalized  weighting  function  used  in  averaging 

I  =  Identity  matrix 

th  =  Gas  phase  mass  production  per  unit  volume 

P  =  Pressure 

Pg  =  Gas  phase  pressure 

Pj  =  Solid  phase  pressure 

aP  =  Phase  stress  difference 

q,  =  Solid  phase  heat  flux 

R  =  Intergranular  stress  (hydrastatic) 

R  =  Irttetgranular  stress  tensor 

r^,  =  Burning  rate  of  the  solid  {^ase  prc^llant 

Sp  =  Interphase  surface  area 

Tj  =  Solid  phase  stress  tetrsor 

Ug  =  Solid  i^iase  velocity 

V  =  Volume 

Vg  =  External  source  velocity 

w  =  Relative  velocity  between  the  phases 
Wg  =  Interfacial  velocity  between  the  media 

X  =  Position  variable 

y  =  Dummy  variable  of  spacial  integration 

e  =  Porosity 
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Eg  =  Gas  phase  porosity 

Eg  =  Initial  loading  Porosity 

Ej  =  Solid  {rfiase  porosity 

=  Solid  i^ase  mass  generation 
p,  =  Solid  phase  density 

Og  =  Bed  mixture  stress 

02*  =  Turbulent  momentum  source  term 
X  =  Dummy  variable  of  temporal  integration 
g  =  Subscript  referring  to  the  gas  phase 

g  =  Sutecripi  referring  to  the  solid  phase 

*  =  Superscript  referring  to  the  fluctuating  portion  of  the  variable 

=  Intetfacial  surface  separating  the  phases 
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Radford  Army  Ammunition  Plant 
ATTN:  SMCAR-QA/HI  LIB 
Radford.  VA  24141-0298 
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1  Commander 

U.S.  Army  Foreign  Science  and 
Technology  Center 
ATTN;  AMXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901-5396 

2  Commandant 

U.S.  Army  Field  Artillery 
Center  and  School 
ATTN:  ATSF-CO-MW.  E.  Dublisky 
ATSF-CN.  P.  Gross 
Ft  Sill,  OK  73503-5600 

1  Commandant 

U.S.  Army  Armor  School 

ATTN:  ATZK-CD-MS.  M.  Falkovitch 

Armor  Agency 

Fort  Knox,  KY  40121-5215 

1  U.S.  Army  European  Research  Office 

ATTN;  Dr.  Roy  E.  Richenbach 

Box  65 

FPO  New  York  09510-1500 

2  Commander 

Naval  Sea  Systems  Command 

ATTN:  SEA62R 
SEA  64 

Washington,  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954-Tech  Library 
Washington,  DC  20360 

4  Commander 

Naval  Research  Laboratory 
ATTN:  Technical  Library 

Code  4410,  K.  Kailasanate 
J.  Boris 
E.  Oran 

Washington.  DC  20375-5000 

1  Office  of  Naval  Research 

ATTN:  Code  473,  R.S.  Miller 
800  N.  Quincy  Street 
Arlington,  VA  22217-9999 

1  Office  of  Naval  Technology 

ATTN:  ONT-213.  D.  Siegel 
800  N.  Quincy  St. 

Arlington.  VA  22217-5000 


No.  of 

Copies  Organization 

4  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  730 
Code  R-13. 

R.  Bernecker 
H.  Sandusky 

Silver  Spring.  MD  20903-5000 

7  Commander 

Naval  Surface  Warfare  Center 
ATTN:  T.C.  Smith 
K.  Rice 

S.  Mitchell 
S.  Peters 
J.  Consaga 
C.  Gotzmer 
Technical  Library 

Indian  Head.  MD  20640-5000 

5  Commander 

Naval  Surface  Warfare  Center 
ATTN:  CodeG30. 

Code  G32, 

Code  G33,  J.L  East 

T.  Doran 

Code  E23  Technical  Library 
Dahlgren,  VA  22448-5000 

5  Commander 

Naval  Air  Warfare  Center 
ATTN:  Code  388.  C.F.  Price 
T.  Boggs 
Code  3895,  T.  Parr 
R.  Derr 

Information  Science  Division 
China  Lake.  CA  93555-6001 

1  Commanding  Officer 

Naval  Underwater  Systems  Center 
ATTN:  Code  5B331.  Technical  Library 
Newport.  Rl  02840 

1  AFOSR/NA 

ATTN:  J.  Tishkoff 

Bolling  AFB,  D.C.  20332-6448 

1  OLAC  PITTSTL 

ATTN:  D.  Shiplett 
Edwards  AFB,  CA  93523-5000 
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3  ALyLSCF 

ATTN:  J.  Levine 
L  Quinn 
T.  Edwards 

Edwards  AFB.  CA  93523-5000 

1  WUMNAA 

ATTN:  B.  Simpson 
Eglin  AFB,  FL  32542-5434 

1  WtyMNME 

Energetic  Materiais  Branch 
2306  Perimeter  Rd. 

STE9 

Eglin  AFB.  FL  32542-5910 

1  WL/MNSH 

ATTN:  R.  Drabczuk 
Eglin  AFB.  FL  32542-5434 

2  NASA  Langley  Research  Center 

ATTN:  M.S.  408.  W.  Scallion 

D.  Witcofski 
Hampton.  VA  23605 

1  Central  Intelligence  Agency 

Office  of  the  Central  References 
Dissemination  Branch 
Rnom  GE-47,  HQS 
Washington,  DC  20502 

1  Central  Intelligence  Agency 

ATTN:  J.  Backofen 
NHB,  Room  5N01 
Washington,  DC  20505 

1  SDIO/TNI 

ATTN:  L.H.  Caveny 
Pentagon 

Washington.  DC  20301-7100 

1  SDIO/DA 
ATTN:  E.  Gerry 
Pentagon 

Washington,  DC  21301-7100 

2  HQ  DNA 
ATTN:  D.  Lewis 

A.  Fahey 

6801  Telegraph  Rd. 

Alexandria.  VA  22310-3398 


No.  of 

Copies  Organization 
1  Director 

Sandia  National  Laboratories 
Energetic  Materials  &  Fluid  Mechanics 
Department,  1512 
ATTN:  M.  Baer 
P.O.  Box  5800 
Albuquerque.  NM  87185 

1  Director 

Sandia  National  Laboratories 
Combustion  Research  Facility 
ATTN:  R.  Carling 
Livermore,  CA  94551-0469 

4  Director 

Lawrence  Livermore  National 
Laboratory 
ATTN:  L-355. 

A.  Buckingham 
G.  Benedetti 
M.  Finger 

L-324,  M.  Constantino 
P.O.  Box  808 

Livermore,  CA  94550-0622 

2  Director 

Los  Alamos  Scientific  Lab 
ATTN:  T3/D.  Butler 

M.  Division/B.  Craig 
P.O.  Box  1663 
Los  Alamos,  NM  87544 

2  Battelle  Columbus  Laboratories 

ATTN;  TACTEC  Library,  J.N.  Huggins 
V.  Levin 
505  King  Avenue 
Columbus,  OH  43201-2693 

1  Battelle  PNL 

ATTN:  Mr.  Mark  Garnich 
P.O.  Box  999 
Richland,  WA  99352 

1  Institute  of  Gas  Technology 

ATTN;  D.  Gidaspow 
3424  S.  State  Street 
Chicago,  IL  60616-3896 

1  Institute  for  Advanced  Technology 

ATTN:  T.M.  Krehne 
The  University  of  Texas  of  Austin 
4030-2  W.  Braker  Lane 
Austin.  TX  78759-5329 
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2  CPIA  -  JHU 

ATTN:  Hary  J.  Hoffman 
T.  Christian 

1 0630  Little  Patuxent  Parkway 
Suite  202 

Columbia,  MD  21044-3200 

1  Brigham  Young  University 

Department  of  Chemical  Engineering 
ATTN:  M.  Beckstead 
Provo.  UT  84601 

1  Jet  Propulsion  Laboratory 

California  Institute  of  Technology 
ATTN:  L.D.  Strand.  MS  125/224 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 

1  California  Institute  of  Technology 

204  Karman  Lab 
Main  Stop  301  -46 
ATTN:  F.E.C.  Culick 
1201  E.  California  Street 
Pasadena,  CA  91109 

3  Georgia  Institute  of  Technology 

School  of  Aerospace  Engineering 
ATTN:  B.T.  Zim 
E.  Price 
W.C.  Strahle 
Atlanta.  GA  30332 

1  Massachusetts  Institute  of  Technology 
Department  of  Mechanical  Engineering 
ATTN:  T.  Toong 

77  Massachusetts  Avenue 
Cambridge.  MA  02139-4307 

2  University  of  Illinois 
Department  of  Mochanical/Industry 

Engineering 
ATTN:  H.  Krier 
R.  Beddini 

144  MEB;  1206  N.  Green  St. 

Urbana,  IL  61801-2978 

1  University  of  Maryland 

ATTN:  Dr.  J.D.  Anderson 
College  Park.  MD  20740 


1  University  of  Massachusetts 

Department  of  Mechanical  Engineering 
ATTN:  K.  Jakus 
Amherst.  MA  01002-0014 

1  University  of  Minnesota 

Department  of  Mechanical  Engineering 
ATTN:  E.  Fletcher 
Minneapolis,  MN  55414-3368 

3  Pennsylvania  State  University 

Department  of  Mechanical  Engineering 
ATTN:  V.  Yang 
K.  Kuo 
C.  Merkle 

University  Park,  PA  16802-7501 

1  Rensselaer  Polytechnic  Institute 

Department  of  Mathematics 
Troy.  NY  12181 

1  Stevens  Institute  of  Technology 

Davidson  Laboratory 
ATTN:  R.  McAlevy  III 
Castle  Point  Station 
Hoboken,  NJ  07030-5907 

1  Rutgers  University 

Department  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  S.  Temkin 
University  Heights  Campus 
New  Brunswick,  NJ  08903 

1  University  of  Southern  California 

Mechanical  Engineering  Department 
ATTN:  0HE200,  M.  Gerstein 
Los  Angeles,  CA  90089-5199 

1  University  of  Utah 

Department  of  Chemical  Engineering 
ATTN:  A.  Baer 

Salt  Lake  City.  UT  84112-1194 

1  Washington  State  University 

Department  of  Mechanical  Engineering 
ATTN:  C.T.  Crowe 
Pullman,  WA  99163-5201 

1  AFELM,  The  Rand  Corporation 

ATTN:  Library  D 
1700  Main  Street 
Santa  Monica,  CA  90401  -3297 
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1  Arrow  Technology  Associates.  Inc. 
ATTN;  W.  Hathaway 

P.O.  Box  4218 

South  Burlington,  VT  05401-0042 

3  AAI  Corporation 
ATTN:  J.  Hebert 

J.  Frankie 
D.  Cleveland 
P.O.  Box  126 

Hunt  Valley.  MD  21030-0126 

2  Alliant  Techsystems,  Inc. 

ATTN:  R.E.  Tompkins 

J.  Kennedy 
7225  Northland  Dr. 

Brooklyn  Park.  MN  55428 

1  AVCO  Everett  Research  Laboratory 
ATTN;  D.  Stickler 

2385  Revere  Beach  Parkway 
Everett.  MA  02149-5936 

1  General  Applied  Sciences  Lab 

ATTN:  J.  Erdos 
77  Raynor  Ave. 

Ronkonkama,  NY  11779-6649 

1  General  Electric  Company 

Tactical  System  Department 
ATTN:  J.  Mandzy 
1 00  Plastics  Ave. 

Pittsfield,  MA  01201-3698 

1  IITRI 

ATTN;  M.J.  Klein 
10  W.  35th  Street 
Chicago.  IL  60616-3799 

4  Hercules,  Inc. 

Radford  Army  Ammunition  Plant 
ATTN:  L.  Gizzi 

D.A.  Worrell 
W.J.  Worrell 
C.  Chandler 

Radford.  VA  24141-0299 

2  Hercules,  Inc. 

Allegheny  Ballistics  Laboratory 
ATTN;  William  B.  Walkup 

Thomas  F.  Farabaugh 
P.O.  Box  210 

Rocket  Center.  WV  26726 
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1  Hercules,  Inc. 

Aerospace 
ATTN;  R.  Cartwright 
100  Howard  Blvd. 

Kenville,  NJ  07847 

1  Hercules.  Inc. 

Hercules  Plaza 
ATTN:  B.M.  Riggleman 
Wilmington,  DE  19894 

1  MBR  Research  Inc. 

ATTN;  Dr.  Moshe  Ben-Reuven 
601  Ewing  St.,  Suite  C-22 
Princeton,  NJ  08540 

1  Olin  Corporation 

Badger  Army  Ammunition  Plant 
ATTN:  F.E.  Wolf 
Baraboo,  Wl  53913 

3  Olin  Ordnance 

ATTN:  E.J.  Kirschke 
A.F.  Gonzalez 
D.W.  Worthington 
P.O.  Box  222 

St.  Marks,  FL  32355-0222 

1  Olin  Ordnance 

ATTN:  H.A.  McEtroy 
10101  9th  Street.  North 
St.  Petersburg.  FL  33716 

1  Paul  Gough  Associates,  Inc. 

ATTN;  P.S.  Gough 
1048  South  St. 

Portsmouth.  NH  03801-5423 

1  Physics  International  Library 
ATTN;  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro,  CA  94577-0599 

2  Princeton  Combustion  Research 

Laboratories,  Inc. 

ATTN:  N.  Mer 

N.A.  Messina 
Princeton  Corporate  Plaza 
11  Deerpark  Dr.,  Bldg  IV,  Suite  119 
Monmouth  Junction,  NJ  08852 
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Qroanigation 

Rockwell  International 
Rocketdyne  Division 
ATTN:  BAGS. 

J.  Flanagan 
J.  Gray 
R.B.  Edelman 
6633  Canoga  Avenue 
Canoga  Park,  CA  91303-2703 

Rockwell  International  Science  Center 
ATT  N:  Dr.  S.  Chakravarthy 
Dr.  S.  Palaniswamy 
1049  Camino  Dos  Rios 
P.O.  Box  1085 
Thousand  Oaks,  CA  91360 

Southwest  Research  Institute 

ATTN:  J.P.  Riegel 

6220  Cutebra  Road 

P.O.  Drawer  28510 

San  Antonio.  TX  78228-0510 

Sverdrup  Technology,  Inc. 

ATTN:  Dr.  John  Deur 
2001  Aerospace  Parkway 
Brook  Park,  OH  44142 

Thiokol  Corporation 
Elkton  Division 
ATTN:  R.  Wilier 
R.  Biddle 
Tech  Library 
P.O.  Box  241 
Elkton.  MD  21921-0241 
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1  Vehtay  Technology.  Inc. 

ATTN:  E.  Fisher 
4845  Millersport  Hwy. 

East  Amherst.  NY  14501-0305 

1  Universal  Prr^ulsion  Company 

ATTN:  H.J.  McSpadden 
25401  North  Central  Ave. 

Phoenix,  AZ  85027-7837 

1  SRI  International 

Prc^ulsion  Sciences  Division 

ATTN:  Tech  Library 
333  Ravenwood  Avenue 
Menlo  Park.  CA  94025-3493 

Aberdeen  Proving  Ground 

1  Cdr,  USACSTA 

ATTN:  STECS-PO/R.  Hendricksen 
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1  Ernst-Mach-lnstitut 

ATTN;  Dr.  R.  Reiser 
Haupstrasse  18 
Weil  am  Rheim 
Germany 

1  Defence  Research  Agency,  Military 

Division 

ATTN;  C.  Woodley 
RARDE  Fort  Halstead 
Sevenoaks,  Kent,  TNi4  7BP 
England 

1  School  of  Mechanical.  Materials,  and 

Civil  Engineering 
ATTN;  Dr.  Bryan  Lawton 
Royal  Military  College  of  Science 
Shrivenham,  Swindon,  Wiltshire. 

SN6  8LA 
England 
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2  instilut  Sami  Louis 

ATTN  Dr.  Marc  Giraud 

Dr  Gunther  Sheets 
Postfach  1260 
7858  Weaii  am  Rhem  1 
Germany 

1  Explosive  Ordnance  Division 

ATTN:  A  Witdegger  Gaissmaier 
Defence  Soence  and  Technology 
Organisation 
P  0  Box  1750 

Salisbury,  South  Australia  5108 

1  Armaments  Division 

ATTN  Dr,  J  Lavigne 
Defence  Research  EstaWishmem 
Vaicanier 

2459.  Pie  XI  Btvd  ,  North 
P  O  Box  8800 

Courceiette,  Quebec  GOiA  i  Ro 
Canada 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undenaJces  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes  Your 
comments/answers  to  the  iiems/qucsiions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ABL-MR-A? _ Date  of  Report  R»rch  1993  _ _ 

2.  Date  Report  Received _ _ _ _ _ _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  imcrtsi  for 

which  the  report  will  be  used.) _ _ _ _ _ _ _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  sourtx,  design  data,  procedure,  source  of 
ideas,  etc. ) _ _ _ _ _ _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man  hours  or  doilars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc'’  If  so.  please  elaborate _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports’’  (Indicate 
changc,s  tc  o.ganization,  technical  content,  format,  etc  ) _ _ _ _ _ _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box.  No. 


City.  Slate.  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Addrc.ss  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


Department  of  the  Army 


OFRCIAL  BUSINESS 


II 

BL'SINESS  REPLY  »L\IL 

nnsT  OASS  femt  m  mai.  m.  m 


Postage  «ill  be  paid  by  addiessee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground.  MD  21005-5066 


